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Abstract. Results of a CCD study of the variability of the cataclysmic variable BG CMi obtained at the Korean 
1.8m telescope in 2002-2005 are presented. The "multi-comparison star" method had been applied for better 
in , accuracy estimates. The linear ephemeris based on 19 mean maxima for 2002-2005 is HJD 2453105.31448(6) + 

0.01057257716(198) (£■ - 764707). The period differs from that predicted by the quadratic ephemeris by Pych et 
al. (1996) leading to a possible cycle miscount. The statistically optimal ephemeris is a fourth-order polynomial, 
as a quadratic or even a cubic ephemeris leads to unaceptably large residuals; Min.HJD=2445020. 28095(28) + 
0.0105729609(57)^; - 1.58(32) • 10"^^ .B^ - 5.81(64) ■ lO^^^S^ + 4.92(41) • lO'^^E^. Thus the rate of the spin-up 
of the white dwarf is decreasing. An alternative explanation is that the spin-up has been stopped during recent 
. years. The deviations between the amplutudes of the spin variability in V and R, as well as between phases are 

' not statistically significant. However, the orbital light curves exhibit distinct difi'erence; the corresponding color 

l/-^ , index shows a nearly sinusoidal shape with a maximum at orbital phase ~ 0.2. The variations of the amplitude of 

■ spin waves shows a short maximum at the phase of the orbital dip. The corrected ephemeris for orbital minima 

is Min.HJD=2448368.7225(12) + 0*? 13474841 (6) • {E - 24849) with a narrow dip occuring 0.07P later. The rate 
of the spin period variation seems to be changed, justifying the necessity of regular observations of intermediate 
polars. 
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^> 1. Introduction disk and a fast flipping of the stream from one active pole 

5^ to another. 

d The object 3A 0729 + 103 has been identified as an inter- 

■ " " ' mediate polar-type cataclysmic variable by McHardy et al. Patterson and Thomas (1993) suggested the double 

(1984). They found night-to-night variations of the mean period of 1827 sec instead of the observed 913 sec. Such a 

brightness with an amplitude up to 0"22 in their sample, model corresponds to nearly equal accreting columns lying 

The orbital period was found to be 3'^14™1 and maximal near the equatorial plane. 

amplitude up to a"15. However, the most striking fea- ^oj-e recently, dc Martino et al. (1995) concluded that 

ture is the 15.2 min periodic variations of variable shape, tj^jg pg^od is the spin one. From a Doppler tomogram 

which may be interpreted as the spin period of a magnetic analysis, HeUier (1999) reported on a "weak spin wave" 

white dwarf. In the General Catalogue of Variable Stars, -with the same period 
the designation of the star is BG CMi. 

The general model for intermediate polars is a red The circular polarization has been detected by Penning 

dwarf filling its Roche lobe, and a white dwarf, the mag- (^986). However, no evidence on its modulation with 

netic field of which is strong enough to disrupt accre- f'"" P^^^^ been detected, contrary to e.g. V405 Aqr 

tion disk completely or at least in its internal parts (cf. (Shakhovskoy and Kolesnikov 1997), which could justify 

Patterson 1994, Warner 1995, HeUier 2001). Norton et al. suggested spin period value. 

(1992) suggested that the 913 sec period is a spin-orbital The 15 min variations show a drastic period de- 
beat period, thus indicating a weak or absent accretion crease, as was originally found by Augusteijn et al. (1991) 
and Singh et al. (1991), and studied in more detail by 
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Table 1. Journal of observations: begin and end of the 
run, the number of observations n, the mean values (m) 



and r.m.s. deviations of the 
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Garlick et al. (1994) even suggested a cubic term in the 
ephemeris. 

The study of rotational evolution of the white dwarf 
needs continuous monitoring of such objects. In the 
present paper we report on further studies of the variabil- 
ity of the star based on the CCD observations obtained in 
the V and R filters at the 1.8m telescope in 2002-2005. 

2. Observations and comparison stars 

The observations have been obtained with a thinned 
SITe 2k CCD camera attached to the 1.8m telescope 
of the Korea Astronomy Observatory and Space Science 
Institute (Bohyunsan) . The instrumental V and R systems 
have been used. To determine instrumental magnitudes of 
stars, the IRAF/DAOPHOT package (Massey & Davis 
1992) has been used. 

Altogether, 1 236 V and 170 R observations have been 
obtained during 9 nights from December 26, 2002 to 
February 3, 2005 (JD 2452635-3405) with an integration 
time 15" in 2002 and 30" in 2004 and 2005. 

The journal of observations is presented in Table 1. 
In an addition to the mean magnitudes (m) in the instru- 
mental system m (see description below), we list the r.m.s. 
deviations am from the mean as an estimate of the char- 
acteristic total amplitude consisting of the variability at a 
few time scales and of the observational noise. 

For further analysis of the individual runs, we mark 
them by the filter letter and five last digits of the integer 
part of the Helicentric Julian Date (HJD) of the beginning 
of the run. 

2.1. Using few comparison stars 

To obtain better accuracy, we have used few comparison 
stars in the field. Their BV magnitudes have been esti- 
mated by Henden and Honeycutt (1995). 
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Fig. 1. The dependence on brightness of the r.m.s. devia- 
tion of stellar magnitudes of the variable and the compari- 
son stars on the corresponding sample mean values. Open 
circles - for the night 53035 with relatively good weather, 
open triangles - for the night 53026 with highly variable 
atmospheric transparency. The curved lines correspond to 
the best fit assuming the Poisson noise for the counts of 
the star and background. Horizontal lines correspond to 
the estimated accuracy of the "artificial" mean weighted 
comparison star. The bottom curve in these two pairs cor- 
respond to better atmospheric conditions and thus smaller 
error estimates. Please note that the brightness is cali- 
brated using results from the night 53035, thus the points 
for the night 53036 with worse weather arc shifted to- 
wards larger apparent magnitude and larger error esti- 
mates. Vertical lines project mean magnitude of BG CMi 
onto corresponding " accuracy-brightness" dependence. 



The procedure of the "artificial" mean weighted star 
has been used, which has been described by Andronov 
and Baklanov (2004) and applied to the comparison stars 
in the vicinity of BG CMi by Kim et al. (2004). The in- 
strumental magnitudes of all comparison stars have been 
checked, and the outstanding values have been excluded 
from the analysis. From the images with all "good" mea- 
surements of the comparison stars, the mean magnitude 
differences have been computed for all comparison stars, 
then the magnitudes of all comparison stars have been 
determined using the magnitude of the "main" compar- 
ison star "05" (Henden & Honeycutt 1995) V=12"^457, 
B-V=0^707. This star is also marked as "Comp 1" by 
Pych et al (1996). Then the weights have been determined 
for each comparison star using iterations as described by 
Kim et al. (2004). Then, for each image, the brightness 
of the comparison star has been determined as the mean 
weighted value from estimates based on the magnitude 
differences between the variable and all comparison stars. 

In Fig. 1, the " accuracy- magnitude" dependence is 
shown for two nights with "good" and "worse" atmo- 
spheric conditions. The corresponding error estimates for 
the "artificial" comparison star are 0^00072 and 0^00114. 
The error estimates of the variable star itself for these 
nights may be estimated from a best non-linear fit as 
^"0067 and ^"0117, respectively. 
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Table 2. Characteristics of the multi-sinusoidal fit of individual runs of observations of BG CMi. The flux is expressed 
in units of 10~^^ erg s~^cm~^A~^. The calibration in V has been made using the magnitude of the comparison star. 
For R, the standard magnitude is unknown and has been arbitrarily set to zero. This makes values of intensity to be 
apparently large. They may be calibrated after further determination of the R magnitude of the comparison star. For 
bad phase distribution caused by an incomplete coverage of the period in the first two runs, the values are highly 
biased and should not be taken for further analysis. Such biased or statistically insignificant values are marked by a 
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1381: ± 2318 
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46: ±49 
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20: ±20 


7.17: ±6.67 


1.39: ±1.29 
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3.527 ±0.042 


0.294 ±0.055 


0.337 ±0.054 


0.273 ±0.055 


0.071:±0.059 


0.050:±0.059 


0.092:±0.056 


V 53036 


3.796 ±0.035 


0.362 ±0.050 


0.343 ±0.049 


0.167 ±0.050 


0.040:±0.049 


0.075:±0.049 


0.062:±0.049 


V 53052 


3.578 ±0.035 


0.442 ±0.051 


0.179 ±0.050 


0.048:±0.051 


0.068:±0.049 


0.054:±0.047 


0.097:±0.048 


V 53385 


3.467 ±0.100 


0.257:±0.158 


0.296:±0.158 


0.296:±0.153 


0.131:±0.152 


0.059:±0.136 


0.011:±0.143 


R 53385 


272816±6747 


31239 ±10432 


17926:±10145 


20531:±9874 


4393:±9968 


6942:±9139 


1100:±9558 


V 53404 


3.494 ±0.086 


0.406 ±0.124 


0.421 ±0.120 


0.072:±0.121 


0.135:±0.121 


0.100:±0.120 


0.075:±0.120 


R 53404 


272413±6526 


31228 ±9288 


32071 ±9066 


7256:±9077 


6246:±9273 


5251:±9189 


3509:±9264 



The observed characteristics of 17 comparison stars are 
listed in Kim ct al. (2004), the corresponding finding chart 
is available at http://uavso.pochta.ru/bgcmi_2.GIF. The 
differences between our instrumental magnitudes Vm and 
the standard ones Vhh (Henden and Honeycutt 1995) are 
well explained by the color reduction formula 

Vin = Vhh - 0'?^0025(16) ± 0.(11)A(1/ - i?)„ (1) 

(Kim et al. 2004). However, for our photometry, we have 
not taken into account the difference between our instru- 
mental VR systems and the standard ones, as there are 
no published estimates of the R magnitudes of the com- 
parison stars. The amplitudes will be expressed in the in- 
strumental systems without any conversion. 

The light curves obtained for 9 V runs and 4 R runs 
are shown in Fig. 2. They exhibit both orbital and spin 
periodicity, as well as cycle-to-cycle variability, which will 
be studied below. 

3. Separation of spin and orbital variability 

3.1. " Running Sine" fit 

The orbital variations show a wave with an asymmet- 
ric minimum, which was interpreted as a grazing eclipse, 
when the accretion structure is partially eclipsed, but the 
white dwarf is not (cf. HcUicr 1999, 2001). The initial 
ephemeris by McHardy et al. (1984) has been revised 
by Augusteijn et al. (1991) and Patterson and Thomas 
(1993). Pych et al. (1993) have added 5 minima tim- 
ings and corrected this revised ephemeris. Patterson and 
Thomas (1993) distingush the orbital dips, which are close 
in phase to the X-ray dips, and the orbital minima approx- 
imately corresponding to the mean between the descend- 
ing and ascending branches of the eclipse at it's middle 
parts. The dip occurs ~ 0.12P after the "mid-eclipse", as 
one may see in their Fig. 6. To study the shape of the or- 
bital curve, they subtracted the best fitting sinusoid with 
the spin period for each night. 



To take into account the variability of the shape of the 
913-sec signal, a local sinusoidal fit 

XRs{t, to, At) = a-r- cos(w(t - Tq) - 27r(po), (2) 

was applied. Here the data at time interval {{to — At), {to + 
At)) arc taken into ac;count for least squares computation 
of the coefficients. The parameters a (local mean value), 
r (semi-amplitude) and (po (phase of maximum in units 
of the short period) arc dependent both on the mid-time 
to of local sub-interval and the filter half-width At. Here 
ui = 27r/P, P is period, and Tq is the initial epoch. To 
follow cycle-to cycle changes, we have adopted the value 
At = 0.5P. 

General expressions for the statistical properties of 
running approximations for arbitrary basic and weight 
functions have been presented by Andronov (1997), with a 
comparative study of "running parabola" , "running sine" 
and wavelet fits for various weight functions (Andronov 
1999). This method had been applied for symbiotic vari- 
ables UV Aur (Chinarova 1998) and V1329 Cyg (Chochol 
et al. 1999), where shorter variations with variable shape 
and nearly constant period are superimposed on (gener- 
ally aperiodic) long-term trends. 

The most recent ephemeris 

HJDMax = 2445020.2800 + 0.010572992 • E (3) 

(Pych et al. 1996) has been used for reference. The long- 
term variability of the period and thus the initial epoch 
discussed in next section is negligible for such a iocai fit, 
the fit obtained for all interval of observations. Thus this 
initial value of the period may be used for local fits. 

As an approximation to the longer-term orbital light 
curve, we have used local values of a{to. At) at times to co- 
inciding with times of real observations. The correspond- 
ing values are plotted versus the orbital phase computed 
using the ephemeris 

HJDMin = 2445020.384 + 0.1347486 • E (4) 
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Fig. 2. The brightness variations of BG CMi during our observations with orbital phase computed using the ephemeris 
(4) by Patterson and Thomas (1993). Left: original data; Right: "running sine" fit Xiis{t, t, At) for times of observations 
and corresponding local mean a{t,t,At). In the fits, some points have been skipped because of large error estimates 
of the smoothing value (> 0™10) near the gaps of the observations. 



(Patterson and Thomas 1993). The r.m.s value of accu- 
racy estimate of a is aa = CK"0147, so total amplitude of 
variations reaches 22aa- 

The corresponding fits for slow variability a{t, t, At) 
and fast+slow variability X/{s(t, t, At) computed for times 
of observations, where these fits are available, are shown 
in the right part of Fig. 2. One may note a significant shift 
of the mid-eclipse as compared to the ephemeris (4). 

Prom our observations, we have determined 3 mo- 
ments of minima using the "running sine" approxima- 
tion, which are listed in Table 4 (runs V53035-V53053). 
Unfortunately, other runs do not cover eclipses completely, 
thus this method was not usable. For those runs (52669- 
52673, 53386, 53404), we have appHed another method 
based on multi-sinusoidal fits, which will be discussed be- 
low. 



3.2. Multi-sinusoidal fits 

De Martino et al. (1995) have applied a multi-sinusoidal fit 
assuming multi-harmonic shape of the orbital light curve 
with frequencies uij = jQ, where Q = 1/P is the fre- 
quency corresponding to the orbital period P, j = l..s. In 
addition, they took into account the 913-scc pulses with 
frequency ujg+i assuming its nearly sinusoidal shape. 

The mathematical model for this two-frequency pro- 
cess is 

s 

I{t) =Io+J2 cos{io,{t - Toj)). (5) 

Here I{t) is intensity (computed from the magnitudes) and 
Aj are semi-amplitudes of variations, and Toj arc initial 
epochs of the maxima of the waves corresponding to fre- 
quencies Uj. We use fluxes instead of magnitudes, to make 
results comparable to those of de Martino et al. (1995). 
According to the calibration (cf. Johnson 1955, Allen 
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Fig. 3. The multi-sinusoidal fit for individual runs of ob- 
servations of BG CMi. The orbital light curve is shown 

as an abridged sum (s = 6), whereas the "orbital+spin" 
variations are fitted by a complete sum (s = 7). The fits 
have been computed for fluxes, following de Martino et al. 
(1995), and converted to stellar magnitudes for compari- 
son with the original data and the "running sine" fits. For 
R, the magnitude is expressed as the "Var-Comp" differ- 
ence. 

1973), the magnitude V=15™00 corresponds to the value 
of the monochromatic flux fx = io-i4.44 ^ 3 . ;^q-15 
erg s~"'^cm~^A~-'^). 

Results of the fit for the value s = 6-1-1 (for orbital and 
spin variability, respectively), which had been adopted by 
de Martino et al. (1995), are listed in Tables 2 and 3. It 
should be noted that the parameter Jq generally differs 
from the sample mean value (/). 

For two short runs 52635 and 52638, the distribution of 
observations in the orbital phase is bad, thus the values of 
the parameters are highly biased and thus have no phys- 
ical meaning, despite the fit itself is being good enough. 
The values deviating from zero less than 3a are not statis- 
tically significant, and marked by the symbol " :" as "bad" 
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Fig. 4. Brightness and color variations for original and in- 
terpolated points (dots) and approximations of the orbital 
and "orbitaH-spin" variability. 

ones. For the first run, all values are thus marked as "bad" , 
despite the amplitude for the spin variability having a rea- 
sonable value. For other nights, the estimates of the mean 
semi-amplitude of the spin variations varied in the range 
(0.38-0.46) -10"^^ erg s'^cm'^A"^. This interval slightly 
cxcecids the error estimate, indicating cycle-to-cycle and 
thus night-to-night variability of this parameter. 

One may note that the ratio Porb/Pspin = 12.74 is far 
from an integer value, thus the orbital light curve differs 
from that of the next cycle by a spin phase shift of a 
quarter of period. 

For all nights, the amplitudes of the waves with the 
orbital period Ri, its first harmonic i?2 and the spin pe- 
riod i?7 are statistically significant. The value R3 slightly 
exceeds 3a only for two nights. For other nights, as well 
as for higher harmonics of the orbital period, the ampli- 
tudes i?3, R4, R5, Re are not statistically significant, and 
generally should not be taken into account. However, we 
present results in Tables 2 and 3 for the same "6-1-1" multi- 
sinusoidal model, as has been proposed by de Martino et 
al. (1995). The nndti-sinusoidal fits for fluxes are shown 
in Fig. 3. The fits show the presence of some variations of 
the orbital light curve, if taking into account such a large 
number of harmonics. 

We may conclude, that, for our VR observations ob- 
tained in 2002-2005, the highest harmonics are not signif- 
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Fig. 5. Two-color diagram for the deviations of original 
and interpolated points from orbital fit (dots). The line 
corresponds to the unitary slope AR = AV. The closed 
lines correspond to the orbital fit (shifted by +0™1 in V 
and — 0™1 in R) and deviations of the "spin+orbital" fit 
from the pure "orbital" one. The orbital curve is shifted by 
+0™1, — 0™1 for illustrating purposes. The cross indicates 
±lcr values for accuracy of the smoothed curves in V and 
R. 

icant, contrary to the B observations of de Martino et al. 
(1995) obtained in 1983-1992. For firm conclusions, new 
multi-color observations of the object are needed. 

The direct comparison of the fluxes between this work 
and that by de Martino et al. (1995) is not possible, as 
they have published their results in another photometric 
system (B). The peak-to-peak amplitude (Table 3) of the 
multi-sinusoidal fit of the orbital variability for 5 runs, 
which completely covers the period, varies from 0™38 to 
0™46 in V. However, for the same nights, the amplitudes 
in V and R are practically the same around 0™45. 

Despite a generally good agreement of the multi- 
sinusoidal fit with the observations, showing both a- 
sinusoidal variability of the orbital light curve and the spin 
pulses superimposed, there are some important differences 
to the "running sine" fits. In the "global" multi-sinusoidal 
model, the amplitude of the spin pulses is constant (in 
the "flux" representation), thus its apparent "magnitude" 
variability with phase is caused by the variability of the 
smoothed phase dependence of the orbital flux, and not 
by the physical variability of individual spin pulses. 

3.3. Night-to-night variability 

The mean brightness variations during 9 nights of our ob- 
servations show only minor variations with a peak-to-peak 
amplitude of 0^10, as one may see from the values of Jq 
in Table 2. This value is smaller than that 0™22 reported 
by McHardy et al. (1984). This may be an effect of a 
small sample. However, if assuming that these variations 



are more pronounced in the ultraviolet, and given that 
McHardy et al. (1984) have used unfiltered observations, 
one may interpret such a difference in amplitudes as a re- 
sult of different photometric systems. Unfortunately, we 
have not found other amplitude estimates in the litera- 
ture, thus the question is still open. 

4. Color variability with the phase of orbital and 
spin period 

4.1. Quasi-Simultaneous VR Observations 

Our observations in V and R have been made using alter- 
natively changing fllters. In this case, direct color measure- 
ments are not possible, as BG CMi shows fast variability. 
To make a set of quasi-simiiltaneous observations, the in- 
terpolating or smoothing technique is to be applied. For 
another cataclysmic variable TT Ari, we have used the 
method of "running parabolae" Andronov (1997) to make 
a set of simultaneous smoothed values in 3 colors UBV 
(Tremko et al. 1996, Andronov et al. 1999). That method 
is effective for stars with flickering-dominated variability. 

For BG CMi, the major amplitude of variations is due 
to orbital and spin variability, thus we have applied an 
interpolating procedure to determine brightness at missing 
times using times of observations in another filter. For the 
local sequence VRVRVRV, we have applied a local cubic 
polynomial fit using two nearby observations before and 
after. For regularly spaced observations, this will lead to 
an interpolating formula Xi = (9(a;i_i -I- Xi-f-i) — {xi-3 + 
Xi^3))/16 assuming that indices z — 2, i -I- 2 correspond to 
the same channel as the point i. The statistical error of this 
smoothing value is ctj = (41/64)^/^(72; « O.SOax, where ax 
is an accuracy estimate of the individual signal value. Such 
an approximation is better than a linear one in our case 
of relatively fast oscillations. For real observations with 
nearly the same time interval between subsequent points, 
the weight coefficients have been computed according to 
a polynomial interpolation. Near the borders, where the 
point z ± 2 is located outside a run of equidistantly spaced 
times, we have applied a parabolic fit using the points i =F 
l,i±l,i±3ora linear interpolation using the observations 
i =F 1) * ± 1- No extrapolation has been made to avoid large 
errors. 

As a result of this procedure (see Andronov and 
Baklanov (2004) for a detailed study of statistical proper- 
ties of such approximations), we have obtained 269 pairs of 
"simultaneous" VR observations, for which the instrumen- 
tal color index V-R has been computed. The data points 
which are present only for one color for a trial time, have 
been excluded. 

As unfortunately the brightness of the comparison 
stars in R is unknown, we cannot convert instrumental 
magnitudes and the corresponding color index to the stan- 
dard system. Thus the magnitudes will be expressed as 
Var-Comp, and the colors will be represented with re- 
spect to the comparison star, i.e. our y-R=(Vvar-^comp)- 

var~Rcomp} • 



R 
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4.2. Brightness and color variations 

The corresponding light and color curves for this set 
of "original+interpolated" points are shown in Fig. 4. 
Despite cyclc-to-cyclc variability of the individual spin 
pulses, their characteristic times are much longer than the 
time resolution of our observations, thus the interpolation 
has not biased the shape of the light curve, as one may sec 
by comparing with the data presented at high-harmonic 
fits for intensities at Fig. 3. Because of the relatively small 
amplitude of orbital variations 0™4), the difference in 
the fits of magnitudes and intensities is compared with the 
error estimates, and may only be important at the grazing 
eclipse and subsequent dip. 

The degree of the trigonometric polynomial for the 
multi-sinusoidal fits have been reduced to s = 2+1, taking 
into account that higher harmonics of the orbital variabil- 
ity are not statistically significant (see Table 2). Despite 
the long runs are presented in Fig. 4 separately, the fits 
have been obtained for a long joint run to decrease the 
night-to-night scatter. Thus the approximation of orbital 
variability is the same for both nights. The spin frequency 
is 94.58 cycles/day or 12.7447 cycles per orbital period, 
thus the phase of the spin variations changes by a quar- 
ter of a period every Port- However, in our observations, 
the interval between the orbital minima was apparently 
UlPorb = 1797.003Pspj„, so we have a rare event of coin- 
cidence of both spin and orbital phases in both nights, so 
the fits for the runs arc the same. 

Under such conditions, all differences between the 
curves are due to physical variability of the star and not to 
the effects of dependence on two independent phases. The 
high-amplitude spin pulsation cycle at iporb ~ 0.5 seen in 
the run 53385 has no counterpart for the run 53404, where 
in the vicinity of this phase an opposite phenomenon (am- 
plitude decrease) have been observed. 

The mean magnitudes with respect to the comparison 
star are mov = 2™612(7) and ttior = 2™830(6), so the 
relative color index is equal to niov — ^or = — 0™218(4) 
and corresponds to a slightly larger temperature than the 
comparison star. This color difference is small enough to 
justify use of the star as the comparison one. 

The orbital light curve has two minima separated by 
a half-period. If the phase zero corresponds to the grazing 
eclipse, the secondary minimum may be interpreted by a 
partial eclipse of the illuminated secondary by the disk 
and/or ellipticity of the disk. The semi-amplitudes of the 
contribution with the orbital period are Riv = 0™116(10) 
RiR = CK" 123(9), with nearly the same values for its har- 
monic: i?2v = (ni9, R2R = <riU{8). 

To determinine effective amplitudes in other types 
of cataclysmic variables with aperiodic light curves, 
Andronov et al. (1999) used the r.m.s. deviation ct* of 
the smoothed values from the mean. For periodic vari- 
ables, one may use an effective semi-amplitude i?* = 
{Ri + i?i)"^/^ = 2-^/^(7*. For our data, such values are 
practically identical for two colors: = CK"166(10) and 
R*R = 0™167(9). Such behaviour of the intermediate po- 



lar BG CMi is intermediate between an increase of ampli- 
tude with decreasing wavelength in non-magnetic nova-like 
variables with superhumps (e.g. MV Lyr, Walker (1954), 
TT Ari, Andronov et al. (1999)) and with eclipses (e.g. 
DW UMa, Ostrova et al. (2005)) or with increasing wave- 
length in magnetic classical polars (e.g. AM Her, Szkody 
& Brownlee (1977)). 

The multi-sinusoidal analysis of the color index shows 
that the amplitude corresponding to the orbital period 
(i?i = 0^"036(6)) is much larger than for other frequencies. 
The star becomes most blue at the orbital phase 0.42 and 
most red at 0.82. 

A similar situation is with the semi-amplitudes of the 
spin variations: R^v = 0"154(9) and R^r = ^"147(8), 
so the color excess of the spin variations AV — R = 
— 2.51g(i?3y/ii3/j) = — 0™5(9) is zero within error esti- 
mates. Nearly the same similarity of the spin amplitudes 
is observed in another intermediate polar FO Aqr but with 
much smaller values R3V = 0^087(6) and R^r = ffi"095(4) 
(Andronov et al. 2005b). 

4.3. Two-color diagram 

The two-color diagram for the deviations of the data from 
the orbital fit is shown in Fig. 5. The line of equal shifts 
AR = aAV satisfactorily represents individual points 
with a = 1, and is in excellent agreement with the si- 
nusoidal component. As there is a very small difference 
between the amplitudes discussed in previous section, the 
slope a seems to be slightly less than unity. However, this 
effect is within the error estimates. Also the phase shift 
between oscillations in V and R is not statistically signif- 
icant, thus making the minor axis of the ellipse vanish. 
There is a scatter of individual points of 0^048, which is 
much larger than the accuracy estimate of the individ- 
ual measurement of the star of such brightness (0^01, see 
Fig.l). Such scatter may correspond to a possible flickering 
with much higher amplitude in V than in R, that is usual 
for accretion. This hypothesis may be supported by a the 
result that the V — R color is a miich more high correlated 
with V than with R. The corresponding correlation coeffi- 
cients are r[V - R,V] = 0.5(5) and r[V -R,R]= 0.20(8), 
respectively. This may be also the result of the contribu- 
tion of "strongly colored" spin cycles, when the color vari- 
ations distinctly deviate from the smooth orbital curve. 
However, based on the present material, it is not possible 
to obtain sufficient statistics on the occurence of "colored" 
and "non-colored" pulses. 

Optical, UV and X-Ray fluxes of intermediate polars 
show highly structered accretion geometry (cf. Kim & 
Beuermann 1995, 1996). Our results can be used to model 
such an accretion flow. 

The difference in the shape of the orbital light curves 
in two colors leads to a change of the trajectory at the 
two-color diagram from a line to an "8-like" figure. Such 
difference in behaviour is the observational basis for future 
theoretic modeling of the orbital variability. 
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(1999) and Paterson and Thomas (1993). 



4.4. Periodogram analysis 

To study differences in light curves in different colors, we 
have applied the periodogram analysis to the deviations 
of the individual data points from the "orbital" O — Corb 
and "orbital+spin" O — Corh+spin fits. For some objects, 
there are peaks present at beat frequencies, which are the 
combinations of the orbital frequency fl and the spin fre- 
quency n (cf. Warner 1986, 1995, Wynn and King 1992, 
Patterson 1994, Hellier 2001). 

The test-function S{f) has been used (sec Andronov 
(1994) for details). The deviations O — Corb show the se- 
ries of strong peaks corresponding to the spin frequency 
and its 19-day aliases (the interval between the runs 53404 
and 53385). The peaks near this frequency practically dis- 
appear for O — Corb+spin , and there is no evidence for the 
spin frequency for the color index for both types of devi- 
ations. However, there are strong peaks near 3fl (but not 
exactly "at"). Also in R there is an apparent peak at ~ 55 
cycles/day. We suggest that these apparent peaks arise 
from instabilities of the light curves. Such quasi-periodic 
oscillations with a time-scale of 15-60 minutes are char- 
acteristic for some cataclysmic variables of similar orbital 
periods (e.g. Tremko et al. 1996, Andronov et al. 1999). 
However, our periodograms show no significant peaks at 
frequences, where the peaks have been found in previ- 
ous years by Patterson and Thomas (1993) and Hellier 
(1999). There is no evidence for significant peaks at the 



beat/multiple frequencies, where such peaks are seen in 
other cataclysmic variables. 

4.5. Scalegram analysis 

The scalegram analysis has been carried out using the 
"running parabola" approximation. The "(t(A<)" scale- 
gram (Andronov 1997) shows a shape characteristic for 
a two-component signal. For such a case, is useful the 
"A(Af)" scalegram (Andronov 2003), which is shown in 
Fig. 7 for 2002-2004, when one-color photometry pro- 
vided better time resolution. The most prominent max- 
imum corresponds to the spin variability with an effec- 
tive period of Pspin,eff = 0^01058 and semi-amplitude 
Tspin.efs = 0^139. This value slightly exceeds the mean 
semi-amplitude of the spin variations obtained from multi- 
sinusoidal fits rms,ef} = -2.51g(((/o + l7)/{Io - It))^^^), 

which is equal to 0^127. This excess argues for more rapid 

variations of the individual cycles and their phase shifts. 
For a pure sinusoidal signal, both methods should give 
closer results. 

The second peak at the "A(Ai)" scalegram corre- 
sponds to the first harmonic of the orbital variability with 
P2,e// = 0^061 and r2,e// = 0^083. The third peak in be- 
tween may correspond to the third harmonic of the orbital 
period, because the shape is not sinusoidal. As expected, 
the effective values for the spin period coincides much bet- 
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ter with the values obtained by other methods than that 
for the orbital variability. This is caused either by the du- 
ration of observations comparable with the orbital period 
or by the non-sinusoidal shape of the orbital variability. 

5. 194-minute orbital period 

5.1. Orbital minima and dips 

The orbital minimum has a specific asymmetric shape at 
its bottom, showing a narrow dip occurring by ^ 0.07P 
later than the mid-eclipse, as one may see from Fig. 6 of 
Patterson and Thomas (1993). To study further period 
behavior, we have used 42 minima published in Table 1 
of their paper. It seems that the one corresponding to 
HJD 2445730.443 seems to be a misprint. Following the 
authors, we have checked the original reference of Singh 
et al. (1991) for the orbital maximum and corrected the 
decimal part of this point to 0.120. Then we have collected 
5 minima obtained by Pych et al. (1996). In this work, we 
have determined 3 moments of centers of minima from the 
function a{to) obtained using the "running sine" fit for the 
nights 53035, 53036, 53052, when the coverage of eclipses 
was good. The moments of dips for these nights have been 
determined as those corresponding to the minumum of 
a(to)- For the nights 53384, 53404, the coverage of eclipses 
is not complete, so we have used the multi-periodic fit 
(s = 6-h 1). 

Additionally, we have analyzed 3 nights of unpublished 
unfiltered CCD photometry with a 120-sec integration 
time, which was obtained by Footc (2003) at the 93cm 
telescope in Utah as a part of activity of the Center for 
Backyard Astrophysics (http://cba.phys.columbia.edu). 
From 3 nights, 4 mean minima have been determined using 
multi-sinusoidal "(2-1-1)" fits, one of them (52673.6927) 
corresponds to the mean light curve for all these nights. 
This "white light" photometry is marked as "W" in the 
text. 

All orbital minima determined in this work are listed 
in the Table 4. 

Our minima show smaller scatter than that from the 
tables by Patterson and Thomas (1993) and Pych et al. 
(1996), where the phases of minima range from —0.13 to 
-1-0.22, because of the inhomogeneity of their collection 
of previously published (or even estimated from the light 
curves) minima. Adding our 11 points to their 47 minima, 
the best fit ephemeris has been obtained 

Min.HJD = 2445020.3894+ 0.13474841 -i; 

± 20 ± 6 (6) 

Min.HJD = 2448368.7225 + 0.13474841 • {E ~ 24849) 

± 12 ± 6 (7) 

The second ephemeris corresponds to the minimum error 
estimate of the initial epoch, which is closest to a sample 
mean of times used. 

From our minima timings, we may estimate phase dif- 
ference between the mid-eclipse and orbital dip to range 
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Fig. 7. The dependence of the test function A(At) on the 
filter half-width Ai. The two highest peaks correspond to 
the spin and the harmonic of the orbital variability. 



from 0.04 to 0.1 IP. The mean phase of our 7 VR minima 
with respect to the ephemeris of Patterson and Thomas 
(1993) is = -0?056 ± 0?007; it reaches 8cr and 

thus is to be taken into account. Similar negative 0-C are 
seen for the "W" observations obtained by Foote (2003). 
According to the new ephemeris (6), the mean phase 
{fmin) = — 0'!'012±0'?007 is much close to zero, but there is 
still some deviation. It may be explained by inhomegeneity 
of the previously published data, where both orbital min- 
ima and dips have been included without separation. The 
dips occur at the mean phase {^Pdip) — -l-0?049 ± 0?013, 
according to the ephemeris (6), so the mean phase differ- 
ence is {(fdip) — i^min) — -l-0?062 ± 0?015. There is also a 
suggestion that the minima in R occur by ^ 0?034±0?007 
later than in V. However, it is based only on two nights 
of two-color observations, when the minima were badly 
covered by the observations, thus it may have resulted 
from statistical fluctuations of observational errors. This 
suggestion has to be checked in future studies. 

The peak-to-peak amplitudes of the orbital variations, 
obtained from the multi-periodic fit, are 0'."43 and 0™41 
for the runs V 53385 and V 53404, respectively. For the 
filter R, these values are practically the same: 0™43 and 
0?M1. 

As the orbital periods of cataclysmic variables un- 
dergo evolutionary decrease (cf. Patterson 1984), we have 
checked the 0-C diagram for the presence of a quadratic 
term in an ephemeris. Unfortunately, the scatter of the 
early data compiled by Patterson and Thomas (1993) is 
rather large (r.m.s. value of 0?069). So, despite much bet- 
ter accuracy of the data in their and our work and that by 
Pych et al. (1996), the obtained value of the Porb deviates 
from zero only by 0.9cr, and thus the hypothesis of de- 
crease of the orbital period may not yet be justified from 
all available data spanning more than two decades. 
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Fig. 8. The dependence on the orbital phase of the local 
mean and semi-amplitude of the spin variability. 

5.2. Phase dependence of the pulse-averaged 
brightness and amplitude 

The characteristics of the "running sine" fits are shown in 
Fig. 8. The most important of them are the local mean 
(pulse-averaged brightness) a{t) and scmi-amplitudc r(t) 
of the spin pulse. The orbital variability is prominent with 
the broad (0.2_P) minimum and the lagged dip. The ampli- 
tude shows night-to-night scatter being relatively constant 
at the phase intervals 0.45-0.50. The phases have been 
computed according to the ephemeris (4) by Patterson 
and Thomas (1993). 

At the phases of orbital minimum, the semi-amplitude 
shows a minimum for all nights. However, the dip is char- 
acterized by an abrupt peak of the scmi-amplitudc. The 
peak lags the dip by ~ 0.02P, that may be partially ex- 
plained by an abrupt ascending branch after the dip caus- 
ing an apparent increase of the amplitude. However, this 
phenomenon may not explain the peak completely, as the 
amplitude really increases, as one may check the original 
and smoothed light curves. 

6. Variability of the spin period 

6.1. Determination of times of extrema 

For the analysis of one-color observations in 2002-2004, we 
have used the method of "asymptotic parabolae" (AP), 



200000 400000 600000 800000 
cycle E 

Fig. 9. 0-C values for times of spin maxima calculated 
with the linear ephemeris (Pych et al. 1996) with a cycle 
numbering corresponding to the best 4-th order polyno- 
mial fit. Fits are shown, which correspond to ephemerides 
published by Patterson and Thomas (1993), Garlick et al. 
(1994) and Pych et al. (1996) for polynomials of order 
2, 3 and 2, respectively. The order of our polynomial fits 
with optimal cycle counting are marked by numbers. Our 
parabolic fit is close to that of Pych et al. (1996) and is 
not shown separately. Because of possible cycle miscount 
after the gap, our data are shown with shifts of integer 
number of periods. 

proposed by Marsakova and Andronov (1996) for timings 
of extrema of signals with cycle-to-cycle variability. The 
tables of 58 maxima and 63 minima have been published 
separately by Andronov et al. (2005a). They have used 
these results to obtain 5 nightly and 2 seasonal "mean" 
maxima timings with much better accuracy. The extrema 
timings show a scatter of phases within a relatively large 
"corridor" of ±0?3, except 3 points deviating from the 
mean value toward even larger values. These 3 points have 
been omitted. The mean deviation of the minima from 
the phase 0.5 corresponding to maxima is only -|-0?003 ± 
0?020 and thus is not statistically significant. This is in a 
contradiction to a qualitative suggestion of more abrupt 
brightness increase than decrease (de Martino et al. 1995). 

For the observations obtained in 2005, the use of al- 
ternating filters had led to half the number of data per 
cycle, and thus it was not possible to apply the same 



Kim, Andronov, Park, Jeon: Orbital and Spin Variability of BG CMi 



11 



Table 3. Characteristics of the multi-sinusoidal fit of in- 
dividual runs of observations of BG CMi (continuation of 
Table 2). 
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Table 4. Moments of orbital minima and dips 
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method, which is effective to study characteristics of in- 
dividual cycles. For determination of the individual ex- 
trema in this case, we have used the multi-sinusoidal fit 
(5) for "Var-Comp" magnitude differences. Because of the 
statistical symmetry of spin pulses, we have used only 
the wave with a spin period without any additional har- 
monics, contrary to fitting the shape of the orbital phase 
curve. These timings for earlier one-color observations are 
in good agreement with the results presented by Andronov 
et al. (2005a), but, for uniformity, we have used results ob- 
tained using the multi-sinusoidal fits. 

The corresponding timings of the "mean" maxima are 
listed in Tabic 3. To obtain more accurate results, wc have 
determined "mean" maximum for all 3 nights of the "W" 
observations, as well as for two long runs in V and R. 

6.2. Models of variability 

The variability of the spin period of BG CMi had been de- 
tected by Singh et al. (1991) and Augusteijn et al. (1991). 
Patterson & Thomas (1993) determined a quadratic 



Table 5. Times of mean maxima from the multi- 
sinusoidal fit of individual runs of observations of BG CMi. 
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Table 6. Seasonal ephemerids for spin maxima of BG CMi 
according to the linear Rt ts = Tq + P ■ {E — Eq). 
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ephemeris based on their own and compiled times of max- 
ima. Garlick et al. (1994) suggested even a third-order 
polynomial to fit the O-C diagram. In more recent study, 
Pych et al. (1996) rejected this suggestion, arguing for a 
sufficience of the second-order polynomial. 

For reference, we have used again the linear part of 
the most recent ephemeris (3) published by Pych et al. 
(1996). The long-term variability of the period and thus 
the initial epoch is negligible for such a Jocai fit, contrary 
to the complete interval. 

In Fig. 9, different models of O-C variations are shown. 
We checked the ephemeris by Pych et al. (1996). Despite 
it satisfactory fits to the data obtained in 2002-03 and 
2005, the spin maxima in 2003 occurred out of phase, so 
such a fit is not acceptable. The ephemeris by Patterson & 
Thomas (1993) leads to a shift of the cycle numbering of 
all our data by unity, thus has the same problem with 
the year 2004. All our data may be well fitted by the 
ephemeris by Garlick et al. (1994) assuming a numbering 
shift of 7 cycles. However, it badly fits the data by Pych 
et al. (1996). Thus no fit from previously published ones 
describes the whole data set. 

After such a long gap of 6 years, the main problem 
is the correct cycle numbering. For this purpose, we have 
determined the best values for separate seasons, using the 
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maxima timings from Table 5, and the preliminary cycle 
numbering from the ephemeris by Pych et al. (1996). They 
are listed in Table 6. The weighted mean value had allowed 
us to obtain a corrected difference of the cycle numbers 
between the initial epochs for the first (2002-03) and last 
(2005) seasons. 

Assuming the new linear ephemeris for all our data 
(Table 6), one has to change the initial ephemeris by Pych 
et al. (1996). Assuming the quadratic model for 0-C, from 
the values at the period E = and a,t E = 764707 (our 
epoch for 2002-2004), we get Q = -2.71(2) • lO^^^^ sig- 
nificantly different from the value published by Pych et 
al. (1996). This shows that the "global" parabolic fit for 
all data is not satisfactory. Assuming an additional cubic 
term Q^E^ with fixed values of other 3 parameters, the 
value of = -1-9.7 • 10~^° had been estimated assuming 
that, for a negative P, the second derivative P is positive. 
This indicates a deceleration of the spin-up of the white 
dwarf. This result is based only on the analysis of period 
values at the beginning and end of known data, and is not 
dependent on possible minor cycle miscounts. 

Our linear ephemeris fits well the last timings of Pych 
et al. (1996), so formally one may apply a composite 
"spline-type" model for the 0-C with parts of spin-up and 
constant periods. It resembles the "asymptotic parabola" 
model (Marsakova and Andronov 1996) for extrema tim- 
ings. Such abrupt period changes are often observed in 
classical eclipsing binary stars (e.g. Kreiner et al. 2001). 
Another approach is to apply high-order polynomials (e.g. 
Kalimeris et al. 1994). 

For such polynomial models, we have chosen a free pa- 
rameter which corresponds to the cycle number Ef of one 
of the minima from our sample. We have chosen the initial 
epoch for the last season 2005. The initial cycle number 
was chosen to correspond to the quadratic ephemeris by 
Pych et al. (1996), i.e. 791717. Then it was modified within 
limits of ±10 cycles. For each trial value of Ef and of the 
degree of model, the following iterations were made. At 
first, this point replaced all our 19 timings, so giving high 
weight for this point. The computed best fit polynomial 
ephemeris has been used to determine the cycle numbers 
of all 19 timings. The numbering of the previous 90 points 
remained the same. Then a final polynomial fit was ob- 
tained. The test function was the r.m.s. deviation of the 
phases from the fit uo-c- 

For the parabolic fit, the numbering remained the 
same, showing unique deep minimum at the ao-c{Ef) 
dependence at Ef = 791717. However, the maximal de- 
viations exceded 0.4P, thus the fit is not acceptable. The 
cubic fit has an optimal solution for Ef — 791714, which 
corresponds to the ephemeris 

Max.HJD = 2445020.2791(3) + 0.0105730179(37)^ 

-5.28(13)10"i^£;2 + 2.10(10)10""^;^ (8) 

However, the 4-th order polynomial ior Ef — 791715 gives 
much better fit to our data, with a scatter smaller by 1.5 
times than that corresponding to the cubic fit. Thus it is 



the statistically optimal fit. The corresponding ephemeris 
is presented in the abstract. However, the "spline- type" 
and polynomial models may also not be ruled out. 

For a final determination of the cycle numbering and 
choice of the appropriate model, one needs to analyse the 
spin maxima obtained between 1996 and 2002. However, 
the results of these observations are not published yet. 

Future new observations will be crucial for the study of 
the rotational evolution of the white dwarf in this system, 
as the ephemerides obtained using different models will 
diverge in the next half-year or more. 

7. Conclusions 

— The ephemeris for the orbital minima has been cor- 
rected. 

— The amplitude of spin variability has an abrupt max- 
imum at the phase of the orbital dip after the mid- 
cclipse. 

— The orbital light curves in V and R differ, with the 
color index having its miminimum (highest tempera- 
ture) at the phase 0?42 and maximum at the phase 
0?82. The statistically optimal degree of the trigono- 
metric polynomial is 2, unlike the previously published 
value 6. 

— The spin variations in both colors have the same ef- 
fective amplitudes and colors, contrary to the non- 
magnetic cataclysic variables with either superhumps 
or total eclipses and to the classical polars. 

— The ephemeris for the spin maxima in 2002-2005 has 
been determined. 

— The rate of the spin- up of the white dwarf has been sig- 
nificantly decreased as compared with previous years. 

— The optimal mathematical model for the spin 
ephemeris corresponds to a fourth-order polynomial. 
However, because of the 6-year gap after the previ- 
ously published observations, other cycle numbering 
over the gap may not be excluded. It may correspond 
either to the cubic ephemeris or to the "spline-like" 
ephemeris with the spin-up replaced by a constant pe- 
riod rotation. 

— The obtained parameters may be used for comparison 
with theoretic models of accretion in magnetic cata- 
clysmic variables. The further monitoring is needed for 
studies of rotational evolution of the white dwarf. 
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